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The equivalence of multipartite quantum mixed states under lo- 
cal unitary transformations is studied. A criterion for the equiv- 
^ ; alence of non-degenerate mixed multipartite quantum states un- 

■ der local unitary transformations is presented. 
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Quantum entangled states are playing fundmental roles in quantum information process- 
ing such as quantum computation, quantum teleportation, dense coding, quantum crypto- 
graphic schemes quantum error correction, entanglement swapping, and remote state prepa- 
ration (RSP) etc.. However the theory of quantum entanglement is still far from being 
satisfied. To quantify the degree of entanglement a number of entanglement measures have 
been proposed for bipartite states. Most of these proposed measures of entanglement involve 
extremizations which are difficult to handle analytically. For multipartite case how to give 
a well defined measure is still under discussion. For general mixed states till now we don't 
even have an operational criterion to verify whether a state is separable or not. As a mat- 
ter of fact, the degree of entanglement of a multipartite quantum system remains invariant 
under local unitary transformations of every subsystems. Therefore the quantum states can 
be classified according to the local unitary transformations. Nevertheless an explicit picture 
of the orbits (geometry and topology) under such transformations is not yet ready. We even 
don't have a general (operational) criterion to verify if two mixed states are equivalent or 
not under local unitary transformations. 

One approach in dealing with the equivalence of quantum states under local unitary 
transformations is to find the complete set of invariants under local unitary transformations. 
Two states are equivalent under local unitary transformations if and only if they have the 
same values of all these invariants. The method developed in [1, 2], in principle, allows 
one to compute all the invariants of local unitary transformations, though in general it is 
not operational. The invariants for general two-qubit and three qubits systems have been 
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studied in [3, 4]. In [5] a complete set of invariants is presented for bipartite generic mixed 
states. In [6] a complete set of invariants under local unitary transformations is presented 
for rank-2 and multiplicity free mixed (bipartite) states. The invariants for pure tripartite 
states are also studied in [7, 8]. 

In [9] another approach is presented in investigating the equivalence problem in terms 
of fixed point subgroup and tensor decomposability of certain matrices. The problem is 
reduced to verify whether a certain matrix is rank one or not. A criterion for the equivalence 
of two non-degenerate mixed bipartite states under local unitary transformations has been 
presented. 

It is rather difficult to deal with the equivalence problem of multipartite states in terms 
of invariants approach. The number of invariants increases quickly when the subsystems 
increase. Nevertheless, we find that the approach in [9] can be easily generalized to multi- 
partite case. In this note we investigate the equivalence of multipartite states under local 
unitary transformations by using the method in [9]. An operational criterion for the equiv- 
alence of non-degenerate mixed multipartite states under local unitary transformations is 
presented. 

We first consider tripartite case. Let Hi (resp. H 2 and H 3 ) be an M (resp. N and 
P)-dimensional complex Hilbert space, with | e^) , i = 1,...,M (resp. \fj), j = 1,...,N and 
\g k ), k — 1, P), as an orthonormal basis. A general pure state on Hi®H 2 ® H 3 is of the 
form 

MNP 

l^ = EEE a ijM) ® \fj) ® \g k ), a ijk e c (1) 
i=\ j=i k=i 

MNP 

with the normalization X/ X/ X/ ai 3 ka *ijk = 1 (* denoting complex conjugation). A tripartite 

%=i j=i k=i 

quantum mixed state on Hi <g> H 2 <E> H 3 is described by a density matrix p which can be 
decomposed according to its eigenvalues and eigenvectors: p = J2i=i P where Aj 

are the eigenvalues and |z/j), i — 1, MNP, the corresponding eigenvectors of the form (1). 

Two density matrices p and p' are said to be equivalent under local unitary transforma- 
tions if there exist unitary operators U\ on H x , U 2 on H 2 and U 3 on H 3 such that 

p' = (Ui ® U 2 <g) U 3 )p{U 1 ® U 2 (g) U 3 y. (2) 

For a Hermitian matrix A on Hi ® H 2 <g> H 3 , the set of commuting matrices B such that 
AB = BA is called the commutant of A, denoted as C(A). Obviously C(A) is a subalgebra. 
We call the set of unitary matrices U such that UA = AU the fixed point subgroup of A, 
denoted as Cu(A), which is a subgroup of the unitary group of all unitary matrices. In 
the following we say that a matrix V on Hi <S> H 2 <S> H 3 is tensor decomposable if it can be 
written as V = V 1 ® V 2 ® V 3 for V x e End{H x ), V 2 E End(H 2 ), V 3 G End(H 3 ). 

If two density matrices p and p' are equivalent under local unitary transformations, they 
must have the same set of eigenvalues Aj, i — 1, n. Let X and Y be the unitary matrices 
that diagonalize p and p' respectively, 

p = XAX\ p' = YAY\ (3) 

where A = diag(Xi, A 2 , Xmnp)- 
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[Lemma 1]. Let G be the fixed point unitary subgroup associated with p. Then p' is 
equivalent to p under local unitary transformations if and only if the coset GXY^ contains 
a unitary tensor decomposable matrix. 

[Proof]. Suppose p' = (U x ® U 2 ® U 3 )p(U 1 ®U 2 ® U 3 )\ then p' = YAY^ = YX^pXYl 
Hence YX^pXY^ = ®U 2 ® U 3 )p(U 1 ® U 2 ® U 3 Y, or 

(C/i ® C/ 2 ® t^^V = P(f/i ® U 2 ® f/a) 1 ^- 

That is, (E/x <g> C/ 2 <g> [/g^YX* e C^p), and Q/(p)XYt = yxtCV(p) contains a unitary 
tensor decomposable element Ui®U 2 ®U 3 . 

Conversely, assume GXY^ contains a tensor decomposable element Ui ® U 2 ® U 3 . We 
have then UXY^ — U± ® U 2 ® U 3 , Up = pU and U is unitary. Therefore 

p' = YAY* = (U 1 ®U 2 ®U 3 yUXAXW(U 1 ®U 2 ®U 3 ) 

= (C/i ®U 2 ® U^UpU^Ur ®U 2 ® U 3 ) = (U t ®U 2 ® C/ 3 )V(^i ® ^2 ® C/ 3 )- 

□ 

Let Z be an a matrix on Hi ® H 2 ® H 3 . If we view Z as a matrix on spaces ifi and 
if 2 <8> #3, it is an M x M block matrix with each block of size NP x NP. Its realigned 
matrix Z\\ 23 is defined by 

Zi| 23 = [vec(Z u ), • • • , vec(Z 1M ), • • • , vec(Z m ), • • • , wec(Z MM )]*, 

Taking Z as a matrix on spaces ifi ® if 2 and H 3 , we have that Z is an MiV x MN block 
matrix with each block of size P x P, and the realigned matrix Z 12 \ 3 is of the form 

Z i2 \ 3 = [vec(Z n ), • • • , wec(Zi MiV ), • • • , vec(Z MN1 ), ■■■ , vec(Z MNMN )}\ 

where for any M x N matrix A with entries a^, vec(A) is defined to be 

vec(A) = [an, • • • , ayy, 021, • • • , a2iv, • • • , «mi, • • • , oma^]*- 

The above operation of realignment could be defined in an alternative way, 

(-Z1I23) im,jknp ijk,mnp j (^12|3) ijmn,kp 

It is shown that a matrix V can be expressed as the tensor product of two matrices V\ and 
V 2 , V = Vi® V 2 , if and only if [10] 

V = veciV^veciy-if . (4) 

Moreover [7], for an MN x MN unitary matrix U, if U is a unitarily decomposable matrix, 
then the rank of U is one, r(U) = 1. Conversely if r(C/) = 1, there exists anMxM matrix 
Ui and aniVxiV matrix U 2) such that [/ = U\®U 2 and 

c/^t = ^ = k -i Im w t = ^ = Hn ^ (5) 

where In (resp. Im) denotes the N x N (resp. M x M) identity matrix, k > 0, and C/ is a 
unitary tensor decomposable matrix. 



3 



[Lemma 2]. For an MNP x MNP unitary matrix U on Hi <8> if 2 <8> #3, £/ is a unitary 
decomposable matrix if and only if the ranks of Ui\ 23 , U\i\z are one, i.e., r{Ui\ 23 ) = r{Ui 2 \ 3 ) = 
1. 

[Proof]. Suppose C/ = Ui®U 2 ® U 3 , then f/i| 2 3 = vec(Ui)vec(U 2 <8> C/ 3 )*, C7 12 |3 = t>ec(lq ® 
U 2 )vec(U 3 Y, obviously r(t/ 1 | 23 ) = r(Ui 2 \ 3 ) = 1- 

Conversely, if r(C/i| 2 3) = 1, there exists an M x M matrix C/i and an NP x AfP matrix 
U 23 , such that [/ = JJ\ ® U 23 . Therefore 



U 



12|3 



/ «11^23 \ 
Ul 2 U 23 

\ ummU 23 J 



(6) 



where uij,j — 1, • • • , M are the entries of JJ\. r(Ui 2 \ 3 ) = 1 if and only if r(U 23 ) = 1, that is, 
there exists an N x N matrix U 2 and a P x P matrix U 3) such that U 23 
is a unitary tensor decomposable matrix. 

Let p and p' be two density matrices with orthonormal unitary matrices X and Y as 
given in (3). Set 



U 2 <S> U 3 , and U 
□ 



V = X 



( 














(7) 



V ^n, ) 

where rij, i = 1,2, ...,r, stands for the geometric multiplicity of the eigenvalue Aj of p, 
J2i n r = MNP, A n . are some unitary x nj complex matrices. The conclusion on bipartite 
case [9] is still valid, namely, p and p' are equivalent under local unitary transformations if 
and only if rank r(Vo) = 1 for some unitary matrices A ni . For the case that p has distinct 
eigenvalues, we have 

[Theorem 1]. If p has distinct eigenvalues, p is equivalent to p' under local unitary 
transformations if and only if 

V = XDY\ (8) 

D = diag(e' t91 , e 1 ® 2 , ...,e ldMNP ), contains a unitary tensor decomposable element for some 
Oi E R. 

[Proof]. Let X and Y be the unitary matrices that diagonalize p and p' respectively, 
p = XAX\ p' = YAY\ where A = diag(Xi, A 2 , Xmnp)- 

As p (resp. p') has distinct eigenvalues, any set Xi (resp. Y"i) of unitary eigenvectors 
corresponding to the eigenvalues of p (resp. p') can then be obtained through the following 
equation: Xi = XU (resp. Yi = YV), where the unitary matrix U (resp. V) has the form: 
U = diag(e iai , e ia2 ,...,e iaMNP ) (resp. V = diag(e il3 \ e if32 , e if3MNP )), then XiD^ = 
XUDiVW = XDY\ where D 1 = diag(e ie '\ e ie '^ e ie Mivp), D = diag(e ie \ e i02 , e l9MNP ), 
9i — oti + 6[ — Pi, i = 1, 2, MNP, hence p is locally unitary equivalent to p' if and only if 
XDY\ D = diag(e ldl , e t02 , ...,e t0MNP ), contains a unitary tensor decomposable element for 
some 6j e M. □ 
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As an example let us consider a density matrix on 2 x 2 x 2, 
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(10) 



p' is in fact a PPT entangled edge state [11]. p' and p have the same eigenvalue set. 
They are not degenerate in the case a ^ 6 ^ c ^ 1 or 2 or 1/2. Calculating the unitary 
matrices X and Y that diagonalize p and p', we have p = XAX^, p' = YAY', where 
A = diag(2,0,l/a,a,l/b,b,l/c,c). Denote D = dmg(a\, d2, d3, 4, ds, d6, aV, dg), we have 
the matrix V defined by (8), 



V = 
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/ (-t*! + d 8 )/2 {d l + d s )/2 (di + d 8 )/2 (-di + d 8 )/2 \ 
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= d 3 = 


— C?4 — — C?5 — dg 


— d 7 = d 8 


= hr(V 1]23 ) 



V 

obviously, when d 1 = —d 2 
Hence p and p' are equivalent under local unitary transformations. 

Now we consider the multipartite case. A general pure state on H 1 
of the form 



Ho 



M N k 

I*m) = a hi2-i M 
k=l ifc=l 



e h ) <8> \ f i2 ) <g> 



\di M )i 



H m is 



(11) 



with a 



iii2---tM a iii2--.iM 



1 J I e *l ) J I /«2 ) J 



\g iM ), i k = 1,2,.. .,N k , k = 1,2,.. .,M, the 



corresponding orthonormal basis of complex Hilbert spaces Hi, H 2 , Hm- Two density 
matrices p and p' are said to be equivalent under local unitary transformations if there exist 
unitary operators Ui on Hi, U 2 on H 2 , ■ ■ ■, and Um on Hm such that p' = (Ui ® U 2 <S> • ■ ■ <S> 



Um)p(Ui ® U 2 <8) • • • <g> Um) ■ For any non-degenerate density matrices we have 

[Theorem 2]. Let p and p' be two non-degenerate density matrices on Hi <S> H 2 <S> ■ ■ .®Hm 
and X and Y the unitary matrices that diagonalize p and p' respectively, p = XAX\ p' = 
YNY\ where A = diag(\i, \ 2 , Xn^-Nm)- P and p' are equivalent under local unitary 
transformations if and only if the NiN 2 - ■ ■ N M x N ± N 2 ■ ■ ■ N M unitary matrix V = XDY\ 
D = diag(e Wl , 



e i62 ,...,e i(>N i N i- N M), satisfies 
r{Vi\ 2 ... M ) = r(Vi2|3-M) 



r(Vi 2 ... 



M-1\M) 



where 



{Vl\ 2 ...M)i 1 i' 1 ,i2---i M i' 2 --- i 'M ('wili2— »M ,»i»2 - "'m ' 
(^ // L2|3-M)jii2ii^,i3--jMi3-i' M = (^)ili2 



(Vi 2 ...M-l\M)i 1 i 2 ---i M _ 1 i' 1 i' 2 ---i' M _ 1 ,i M i' M (^) ni2---«M I *i* 2 "- i M ' 

We have studied the equivalence of multipartite quantum mixed states under local uni- 
tary transformations in terms of analysis of fixed point subgroup and tensor decomposability 
of certain matrices. A criterion for the equivalence of non-degenerate mixed multipartite 
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quantum mixed states under local unitary transformations has been presented. In fact this 
approach works for general multipartite mixed states. But then in stead of the rank r(V), 
one has to verify if r{V ) could be one for all possible matrices A ni , which is again compli- 
cated. The problem is dramatically simplified when the degeneracy of the related density 
matrices is reduced. In particular, for the non-degenerate case, two density matrices are 
easily verified whether they are equivalent or not under local unitary transformation. 
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